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I. INTRODUCTION 



Systems which are created in ultra-relativistic A-l-A collisions at SPS and RHIC are initially very dense and display 
collective behavior. The quasi-macroscopic nature of those systems allows to study new states of matter at the phase 
, boundary between a hadron gas (HG) and the quark-gluon plasma (QGP). Among the different methods to study it, 
an important role belongs to precise measurement of hadronic observables. One- and multi-particle hadron spectra, 
however, do not carry direct information about the phase transition since one should expect strong interactions in the 
hadron resonance gas after hadronization. Usually, different theoretical approaches are used to diagnose the matter 
properties at the transition point. A rather fruitful concept is that of chemical freeze-out. It implies, first, that the 
initial hadronic state is in chemical equilibrium [^]-^. Second, the chemical concentrations of HG do not change 
during the evolution. The latter assumption is justified by the observation that the rate of expansion in HG is larger 
than the rate of inelastic reactions and less than that of elastic collisions [|| . 
^ ' It gives the possibility to use the approach of local thermal equilibrium and hydrodynamic expansion of HG at 
(— I , approximately frozen chemical composition. The initial conditions for the hadronic stage of the collective expansion 
are determined by matter evolution in the pre-hadronic phase and also by particle interactions and hadronization 
conditions. 

If chemical freeze-out takes place, it is defined by a temperature and chemical potentials conjugated to conserved 
^ ' quantum numbers. To extract those thermodynamic parameters from total particle number ratios one has to assume 
uniform temperature, baryon and strangeness chemical potentials as well as a unique hadronization hypersurface 
(t'^^{x) for all particles. The result is a common effective volume for all hadron species, that completely absorbs effects 
of flow and a form of the hypersurface and approximately cancels in a description of particle ratios pJq] . 

In detail, the situation may be somewhat more complicated. While baryon stopping seems to be fairly complete 
at AGS energies (see also Ref. [|j where possible deviations from the full stopping are analyzed quantitatively), 
there are clear indications for an onset of transparency at top SPS energy from the observed fairly flat rapidity 
distribution of pions in Pb-|-Pb central collisions [P,p^ and, in particular, from the difference in rapidity distributions 
between protons and antiprotons, which cannot be explained by longitudinal flow developing after full stopping. As 
a consequence, one could expect particle number ratios taken in mid-rapidity, especially involving nucleons, to differ 
somewhat from those integrated over 47r. On the other hand, a consistent hydrodynamic analysis then must be based 
on an analysis of data in a narrow mid-rapidity region, which is done in this paper. We also note that, in a boost- 
invariant scenario as observed at RHIC, chemical freeze-out conditions are the same for any individual rapidity slices 
taken within a wide enough rapidity interval. 

The object of this paper is to reconstruct the hadronization stage within the hydrodynamic approach. We assume 
that hadronic matter at that stage has locally a maximal entropy. Effects of hadron-hadron interaction (excluded 
volumes ||ll|, mean field ||l^) as well as specific conditions of hadronization (e.g., possible conservation of number of 
quarks during the hadronization |l3] ) are rather complicated and model dependent. So it is naturally to use minimal 
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number of parameters and take all that into account in the simplest way, i.e., as a unique factor changing the particle 
numbers. In other words, we describe the system at the stage of hadronization as a mixture of ideal gases with 
additional common chemical potential Hch- Afterwards the system expands hydrodynamically with frozen chemical 
composition, cools down and finally decays at some thermal freeze-out hypersurface a*^{x) which could be different 
for different particle species [p^[l^]. Then long-lived particles stream freely into detectors, resonances decay and also 
contribute to spectra of observed particles. We reconstruct a state of the system at the earlier hadronization stage, 
find energy and particle number densities and also estimate roughly the effects of interaction at the hadronization 
stage using an approach described in the following. 

From the analysis of slopes of transverse momentum spectra and of interferometry radii for different particle species 
one can restore the space-time extensions, collective flow parameters and temperature at the thermal freeze-out stage, 
when the thermal system decouples and particle momentum distributions change no longer. The method of joint 
analysis of the spectra and correlations was used in Refs. |p|,p7|-p^ to reconstruct the thermal freeze-out stage. Note 
that, to evaluate the detailed behavior of the (pion) spectrum in the whole transverse momentum region, one needs in 
detail knowledge of the unobserved resonance yields. They can be estimated from the temperature and baryochemical 

§'ential at chemical freeze-out. An attempt to determine the resonance yields in that way was undertaken in Rcf. 
. There however, a trajectory assuming isentropic expansion in full chemical equilibrium was employed following 
Ref. [p^ , in contradistinction to the assumption of chemical freeze-out. In this paper we make a new analysis of 
the thermal freeze-out stage within the hydrodynamic approach for lead-lead collisions based on recent experimental 
data on one particle spectra and HBT correlations for different hadron species. We estimate the basic parameters 
at thermal freeze-out from an analysis of the slopes of particle spectra and of the correlation functions at regions of 
high enough transverse momentum to minimize the influence of resonance decays. We use also the previous results 
from Refs. to determine the influence of uncertainties in reconstruction of the thermal freeze-out stage on the 

reconstruction of the hadronic gas parameters at chemical freeze-out. 

Recent evidence p]-p|,pO| support strongly the early interpretation that chemical freeze-out at SPS and higher 
energies takes place at or near to the phase boundary (see also the discussion in [Q). We employ a similar ap- 
proach here with emphasis on the analysis of particle ratios near mid-rapidity. To complete the reconstruction of 
the hadronization stage and find its space-time extension, we solve the hydrodynamic equations with the Cauchy 
conditions that are found at thermal freeze-out, and consider the solution backwards in time until the chemical freeze- 
out temperature is reached. The method leads to unique results due to the scale invariance of the hydrodynamic 
equations for a mixture of ideal Boltzmann gases with reference to scale-transformations of pressure and densities of 
energy and particles at a fixed temperature. For such a system the variation of initial densities by a unique factor 
results in the same scale transformation of the solutions. Therefore, the temperature as function of time does not 
depend on the common chemical potential, which defines absolute values of densities: the temperature depends only 
on concentrations of different particles species. The latter do not change during the hydrodynamic evolution and 
are determined by the baryonic and strangeness chemical potentials and temperature just after hadronization. Fi- 
nally, from reconstructed space-time extension of the chemical freeze-out and a particle numbers in mid-rapidity we 
determine the densities at hadronization and the common chemical potential. 

An alternative way to determine hadron gas parameters is used in Ref. |14|. There, a hydrodynamic approach is 
followed from the initial condition of an equilibrated QGP through the phase transition and hadronization. Subsequent 
to hadronization the system is assumed, in fact, to evolve as a mixture of mutually non-interacting gases with fixed 
chemical composition. In the present approach we, instead, solve the hydrodynamic equations for a complete system 
under the explicit condition of particle numbers conservation. The distinctive feature of our method is the analysis 
backwards in time. 

One more important point which we would like to comment here concerns the question of applicability of the 
hydrodynamic approach to post hadronization stage in A-l-A collisions. In recent papers it is stated that the 

hydrodynamic approach is unable to describe the stage of hadron phase expansion in heavy ion collisions. Note, 
however, that the criticism of Refs. jl^,^ as to the hydrodynamics is based on calculations of chemically equilibrium 
evolution of the mixture of hadron gases. This evolution differs essentially from the evolution of chemically frozen gas 
p^ , p5t . In very recent work the first detailed comparisons of RQMD calculations with chemically frozen hydro- 
dynamics were done. There it was demonstrated that if chemical freeze-out is incorporated into the hydrodynamics, 
then the final spectra and fireball lifetimes are insensitive to the temperature at which the switch from hydrodynam- 
ics to cascade is made. It implies that the transformation of heat energy into collective fiow for a chemically frozen 
hydrodynamic evolution results in spectra which correspond approximately to the microscopic cascade calculations. 
Moreover, the initial conditions of hydro-expansion in Refs. [y_6[23| correspond to chemically equilibrium mixture of 
Boltzmann hadronic gases and does not include possible specific effects of hadronization process, excluded volumes 
and collective interaction. The estimation of these effects, which may reducing the particle number, is one of the 
important aims of this paper. Therefore, the criticism in Refs as to applicability of hydrodynamics is based 

on the hydrodynamic models which differ essentially from what is used in our paper. The opposite to Refs. 23 
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conclusion as to applicability of hydrodynamics to the post hadronization stage in A+A collisions, was argued recently 
in Refs. ||2^ ]. There it was shown, in particular, that the hydrodynamics describes good the single particle spectra 
and elliptic flows also in RHIC experiments. 

Nevertheless, the cascade models of post hadronization stage are considered now as favorable alternative 

to the hydrodynamic approach. The picture of evolution and particle emission are rather different in hydrodynamic 
approach and cascade models (RQMD, UrQMD, etc.). While in standard hydrodynamic picture the particle spectra 
are evaluated usually according to the Cooper- Frye prescription (using the freeze-out hypersurfaces), the cas- 
cade calculations does not demonstrate the sharp region of particle emission p^ . It does not mean, however, that 
hydrodynamic description is incorrect: the absence of principal contradictions between the both methods of spectra 
calculations was demonstrated quite recently in Ref. based on exact solution of Boltzmann equation. The method 
of escape probabilities, proposed there, gives possibility to make bridge between hydrodynamic and transport models 
for rarefied systems. As to rather dense systems, similar to those that form just after hadronization, where the time 
of collision is similar to the time between collisions and collective effects, such as multiparticle collisions and modifi- 
cation of hadron properties, may take place, the oversimplified version of transport approach, like UrQMD, cannot 
be applied correctly while hydrodynamic approach is more suitable. Note also that at CERN SPS and RHIC energies 



the cascade models fail actually to describe the interferometry radii of the system (see, e.g., [29 23[|). 

In Sec. II we discuss the possibility to extract the thermodynamic parameters at chemical freeze-out from the 
particle number ratios when those parameters are non-uniformly distributed in rapidity. In Sec. Ill we consider the 
particle number ratios in a hydrodynamic approach to A+A collisions. The optimization procedure is applied to 
fit recent experimental data for Pb-f Pb (CERN, SPS) collisions and find the temperature, baryonic and strangeness 
chemical potentials at chemical freeze-out. In Sec. IV we analyze the property of hydrodynamic expansion in a central 
slice of space-time rapidity for a chemically frozen mixture of ideal Boltzmann gases. Approximate hydrodynamic 
solutions for non-relativistic transverse fiows are found. In Sec.V the Cauchy conditions at thermal freeze-out, such 
as proper time, transverse Gaussian radius parameter, temperature and transverse velocities are defined from the 
analysis of spectra and HBT correlations. In Sec. VI the corresponding " inverse in time" hydrodynamic solution is 
analyzed and the space-time characteristics of the matter at the hadronization stage are reconstructed. The later are 
applied to restore the energy and particle number densities at that stage. The conclusions and a discussion of the 
results are presented in Sec. VII. 



II. THE PARTICLE NUMBERS AT MID-RAPIDITY FOR EXPANDING HG 

We assume that just after hadronization, the system can be described at some hypersurface a^^{x) as a locally 
and chemically equilibrated expanding hadron resonance gas. The effects of interaction and possible peculiarities of 
the hadronizaton process are taken into account by a common factor changing particle numbers in the ideal gas; the 
factor is associated with an additional common chemical potential ^ch- In concrete evaluations, as was performed 
by many authors, instead of considering a finite chemical and thermal freeze-out hypersurfaces (t(x) in Minkovsky 
space, we use the common profile factor p{x) at proper time r(r) « const to guarantee the finiteness of the system 
(see the discussion on the validity of such an approximation in ) . Then the total multiplicity of a particle species 
i at chemical freeze-out is calculated by 

where f{p^iU^^(x)/Tch{x),^i^ch{x)/Tch{x)) is the Bose-Einstein or Fermi-Dirac distribution function in our approach , 
and ^ii^ch = Bi^B + Si^s + A*c/i- The additional chemical potential /ich is unique for all particles species. The effect 
of isotopic chemical potential has been investigated in and found to be small. It is neglected in this paper. Taking 
into account the invariance of the distribution function under Lorentz transformation, one can carry out integrations 
over momentum variables and finally get 

N.,^ j daf{x)u^{x)nt{x), (x) = nf (x)p(x) (2) 

where in Boltzmann approximation for the distribution function the thermodynamic densities ri1^{x) have the form 

nt{x) = ^^;^T,h{x)mi eMt^rMx)/Tch{x))K2{mjTM), (3) 

and Kn{u) = i rfz exp[— ixcoshz -I- nz\ , with Rew > 0, is the modified Bessel function of order n (n = 0, 1, ...). 
One can rewrite (^ in the form 
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= {nf)V'"' (4) 

where the volume of system is expressed through the integral of the hydrodynamic 4 -velocity and the 

profile factor p{x) over the hypersurface a'^^: 



V^^ = j dcjf{x)u^^{x)p{x). (5) 

We further define the density averaged over the chemical freeze-out hypersurface by 

. f dafuf'ix^nfix^pix) 
J da^"-ut^(x)p(x) 

If temperature and chemical potentials are uniform over the chemical freeze-out hypersurface, then (n^'') = 
nf'' [Teh, Pi,ch) and the standard analysis of particle number ratios in 47r geometry can be performed to find the 
temperature and chemical potentials. However, if T and p are not uniformly distributed over the hydrodynamic 
tube at the hypersurface a'^^, then the similar analysis implying the substitution Ijii) rif^{{Tch), (pi.ch)) with some 
average temperature and chemical potentials, could result in significant deviations for the values one is interested 
in pO| . To reduce the corresponding uncertainties, the analysis of particle ratios in a narrow rapidity interval near 
mid-rapidity should be used as we will discuss below. 

To analyze non-uniformly distributed systems, let us start from the well known observation that, at the latest 
stage of a hydrodynamic evolution, the longitudinal velocity distribution corresponds to the asymptotic quasi-inertial 
regime, vl = z/t. This property, which is the ansatz in the Bjorken model jS^] for each moment of time, is realized at 
the freeze-out stage even in the Landau model of complete stopping Consequently, we assume in the following 

the boost-invariance for longitudinal hydrodynamic velocities at the last (hadronic) stage of the evolution, vl = z/t. 
This then allows us to parameterize t and z at ct by the form: t = T(r, 77) cosh(77), z = r(r, ry) sinh(r7), where r is 
absolute value of the transverse coordinate r and i] is the longitudinal fluid rapidity, 77 — tanh^^ v^. The parameter 
T(r, rj) generalizes the Bjorken proper time parameter t. Taking into account also a transverse velocity component 
(f(r, r) in the longitudinally co-moving system, LCMS), we obtain for the 4-velocity 

(r, 77) = 7 (cosh ry, v cos 0, v sin 0, sinh 77) , (7) 

where 7 = (1 — w^)^^/^. The element of the hypersurface (t{x) takes the form 

, 1 (ir dr dr 1 dr . , \ /o\ 

dan = T(r, rj)drjdrxdry — — smhry + cosh?y, — - — , — - — . — coshry — smhry . (8) 

\ T dry drx dry ' t drj J 

Let us, first, introduce for an expanding HG the effective volume V^ff attributed to one unit of rapidity. To 
determine Vef / one can write the particle number distribution differential in fluid rapidity ry for particle of species " i" 
at chemical freeze-out: 

dN^^ I d^p daf^^^ / p,n^(r, ry) M,,,(ry) \ ^^^^ 
dr] J p° drj \ Tchiv) ' Tchirj) J 

Here we assume that the chemical potentials as well as the temperature are constants in transverse direction across the 
freeze-out hypersurface. Taking into account the invariance of the distribution function under Lorentz transformation, 
one can carry out integrations over momentum variables and finally gets 

^-Veffiv)nf{v) (10) 



and 



Veff{v)= I -^u^{r,v)p{r,r]). (11) 



^The asymptotic quasi-inertial property one can see also in the analytic solutions for elliptic flow in non-relativistic approxi- 
mation [K3|. 
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The particle number distribution over momentum rapidity y is 



- ( dal'^rf ( ^ ) Pir, V). (12) 



dy J p°dy V Tchiv) ' Tchiv) 

One can rewrite Eq. ( p^ ) in the form 

dN, dN,, 



\n=ya + S^{y))^K7fiy)il + SM)nriy) (13) 



where 



dNj _ dth I 

From Eq. one can conclude that, for a nonuniform distribution of thermodynamic parameters along the beam 

axis, the analysis of particle number ratios at mid-rapidity, y « can be based on ratios /rij if correction factors 
Si(0) are small. One can include corresponding corrections in the non-unique "effective volume" y^f^fi = ^e'//(l + '^i) 
which is different for different particle species. 

Let us estimate the correction factors. The saddle-point approximation to the integrals over fluid rapidity rj in 
can be applied when the parameter mi/T{ri) 3> 1 at the saddle point rj ^ y. Then, at mid-rapidity, the area 
around t] = y = gives the main contribution to the integral over 77 in ( [l^ ) and one can expand the slowly varying 
functions into a Tailor series around the saddle-point 77 = 0. For the sake of simplicity, we neglect a variation of 
T{ri) in space rapidity in comparison with the variation of ^f^{rj), assuming that all non- uniformity originated from 
incomplete baryonic stopping. This implies an increase of the baryonic chemical potential with the absolute value of 
rapidity. Assuming that the transverse expansion is non-relativistic, with a freeze-out that occurs at approximately 
equal proper time r, and that the profile p{x) can be factorized, p{x) — pT{r)pL{vi), we estimate the correction factor 
at mid-rapidity: 

^ ^ 1 (PL exp{pi^ch/Tch)Tch)" \ri=0 Tch ^-j^g^ 

2 pL{0)GyiY>{pt,ch/Tch)Tch rui 

All thermodynamic parameters here are taken at mid-rapidity. One can see that the factor 5i is rather small for 
baryons, because the ratio Tch/fni is small, as well as it is small for pions, because their baryonic number are equal 
to zero. Therefore, as follows from Eq.(^^ for 5i the analysis of particle number ratios at a narrow mid-rapidity 
interval can be based on ratios These values are the same as for Boltzmann gas in a static volume (5||^. 

It is worth to note that the second derivative in formula ( p^ ) is positive for incomplete baryonic stopping and, so, 
5i > 0. As we will discuss later, the consequences of it are non-trivial and rather important for pion production via 
decay of baryonic resonances after the thermal freeze-out occurs. 



III. THE ANALYSIS OF PARTICLE NUMBER RATIOS 



In ultrarelativistic A-f A collisions a hadron resonance gas is assumed to appear due to the hadronization process. 
Subsequent, it expands and reaches the final stage of the evolution - the thermal freeze-out. According to the 
hypotheses of chemical freeze-out, the total number of each particle species is preserved during the evolution. After 
the hydrodynamic tube decays and produces final particles, one has to take into account also particles that come from 
resonance decays. In our calculations the resonance mass spectrum extends over all mesons and baryons with masses 
below 2 GeV. Decay cascades are also included in the analysis. The value of the strangeness chemical potential ps 
is fixed by baryonic chemical potential and temperature from the condition of local strangeness conservation. We 
optimize the particle number ratios at mid-rapidity, using the corresponding values for Pb -I- Pb collisions at 158 A 
GeV/c from WA97 (Table 1) and NA49 (Table 2) experiments, excluding from optimization procedure pions as the 
lightest particles, the later circumstance can result in serious corrections, besides of Bose-Einstein enhancement effect 
in the distribution function, which we will discuss below. Then we describe the data in the model of a chemically 
equilibrated ideal hadron gas (HG). The feeding in our calculations is tuned according to the experimental conditions 
as described below. 

In the geometry of the WA97 experiment the feed-down from weak cascade decays is expected to be of minor 
importance pi], and so the A (A) have not been corrected for feed-down from S weak decays, and S have not been 
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corrected for feed-down from fl. But feeding of A (A) from electromagnetic decays of T,^ (S ) is included. We did not 

include (according to the NA49 experimental procedure JlO|) feeding of p and p from A, A and from . In the 

NA49 experiment A and A yields include the contributions of and if as well as feed-down from weak S decays . 
The feeding of 7r~ from kaons and from A is not included We included feeding from (j) decay into kaons (NA49 
and WA97) and pions (NA49). For technical reasons in the model calculations we used the following substitution: 
/tt~ {K^ + _ftr+)/27r~, a good approximation for Pb-)-Pb collisions Ba]. The results for the best fit are shown 
in Tables 1 and 2 for Boltzmann statistic (BS) and quantum statistic (QSjcalculations. The obtained temperatures 
Tch ~ 164 MeV and baryochemical potential fiB ~ 224 MeV differ slightly from the corresponding values of Ref. 
Tch ~ 168 MeV, /is « 266 MeV, that were gotten from analysis of particle number ratios in different limited rapidity 
intervals as well as in full 47r-geometry. 

TABLE 1. Experimental particle number ratios (WA97) in Pb-I-Pb SPS collisions , compared to model calculation 
for Boltzmann statistic (BS) and quantum statistic (QS) 



Tch 




163.37 MeV 


163.91 MeV 


PB 




223.6 MeV 


224.0 MeV 






55.2 MeV 


55.1 MeV 


Ratios 


Data 


BS 


QS 


A/A 


0.133 ± 0.007 [ 


34 1 


0.133 


0.133 




0.249 ±0.019 [ 


34 1 


0.250 


0.249 


S-/A 


0.100 ±0.004 [ 


34 1 


0.098 


0.099 


S+/A 


0.188 ±0.016 [ 


34 1 


0.186 


0.186 




0.182 ± 0.022 [ 


34 1 


0.139 


0.140 




0.281 ± 0.053 [ 


34 1 


0.274 


0.274 




0.383 ±0.081 [ 


34 1 


0.492 


0.489 




0.6706 |36[ 


0.6212 


0.6226 



TABLE 2. Experimental particle number ratios (NA49) in Pb+Pb SPS collisions , compared to model calculation 
for Boltzmann statistic (BS) and quantum statistic (QS). 



Tch 




163.37 MeV 


163.91 MeV 


Mb 




223.6 MeV 


224.0 MeV 


MS 




55.2 MeV 


55.1 MeV 


Ratios 


Data 


BS 


QS 


[p - p)/7r" 


0.184 [3 




0.245 


0.239 


{p~p)lK- 


1.776 [3 


35 




1.847 


1.843 


{p-p)/K+ 


1.028 [3 


35 




1.090 


1.086 


K+/K- 


1.727 [] 


5| 


1.695 


1.697 


p/p 


0.115 [3 




0.066 


0.067 


(l>/{p-p) 


0.0802 


|35 37] 


0.0937 


0.0935 



The deviation between calculated particle ratios and experimental data are within the error of the data of WA97, 
excepts for the ratios which contain SI. As to the p/p ratio of NA49, one has take into account that the corresponding 
data p/p « 0.11 are rather preliminary and differ much from published NA44 data where the value is twice less [Q. It 
is interesting to note, that in rapidity interval 2.0 < y < 2.5 shifted by 0.65 units of rapidity from central mid-rapidity 
point, the correspondent value measured by NA49 is 0.06 ± 0.01. 

As one can see from Table 2 the most significant deviation between model and data is for the {p — p) /tt^ ratio where 
the calculations predict a value about 30 percent larger then observed. Considering typical sizes of systematic errors 
in the data of at least 10 percent one should, however, not overestimate the statistical relevance of the discrepancy. 
In part, it maybe due to the reduction, compared to the results of Ref. |^], in ^ib since there are less baryons at mid- 
rapidity then are expected in scenario of full baryon stopping. Because about 1/2 of the pions result from baryonic 
resonance decays this also leads to a reduction of pions in the model and, thereby, over-prediction of the (p — p)/iT^ 
ratio. 

The above effects of the non-uniformity in rapidity of the baryochemical potential can be estimated numerically. 
First, note that our result for [p — p)/7r^ ratio implies the same effective volume for pions as for other heavier 
hadrons. In fact, the pion volume is larger due to pion contribution from resonance decays. To explain that, let us, 
first, note that because of large resonance masses, comparing with the thermal freeze-out temperature, the momentum 
rapidity y of a resonance coincides with the rapidity r] of the corresponding decaying fluid element and, therefore, the 
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distribution of resonances in y follows their distribution in space-time rapidity r). The latter is determined mainly 
by the baryochemical potential that has a shallow minimum at = 0. When resonances decay into the lightest 
particles, pions, there is a rather wide kinematic "window" Ay in rapidity for emitted pions. Then resonances with 
y beyond mid-rapidity may affect the multiplicities of pions at mid-rapidity. It is similar to a diffusion process: 
there is transport of pions from baryon rich peripheral y- regions to the central region. One can roughly estimate 
the effect of pion "diffusion" towards mid- rapidity using an effective temperature Tdec for those pions resulting 
from baryon decays. Taking into account that m-n/Tdec ~ 1, one can see from Eqs.(13) and (15) with substitutions 
PL (ixp{fJ'i,ch/Tch)Tch PL^^Tth, Tch/rrii Tdec/m-TT, that a correction, increasing the effective volume for pions 
produced by baryonic resonances, could appear. Our numerical estimate is based on the baryonic density profile 
for AGS energies taken from Ref. [||, gives a contribution of about 10 % for pion transport to mid-rapidity region. 
Whether there are other mechanisms increasing, in particular, the number of pions at mid-rapidity, such as decay of 
a potentially formed disoriented chiral condensate, cannot be decided by the present analysis. 

To provide the hydrodynamic calculations and estimate the energy density one may, phenomenologically, "fix" 
{p — p)/tt^ ratio by introducing an additional chemical potential for direct pions at chemical freeze-out of /Itt = 89 
MeV. Note, however, that the apparent anomaly in (j> — p)/tt~ ratio has only little influence on the hydrodynamic 
behavior of chemically frozen HG considered in the next section. 



IV. HYDRODYNAMIC EVOLUTION OF CHEMICALLY FROZEN HG 



We describe hadron matter evolution at mid-rapidity after hadronization by the relativistic hydrodynamics equations 
for an ideal fluid: 

a^T^" = 0, (16) 
d^u'^n, = 0, (17) 

where T'"' = (e +p)u^u'^ — g^'^P is the energy-momentum tensor of the Boltzmann-like gas (see Sec. II), e = is 
energy density and p — ^Pi is pressure. The energy densities are proportional to the profile factor p{x) similar to 
density of particle species "i" (^: £i{x) — p{x)ei{x). Here £i(x) are thermodynamic energy densities in the Boltzmann 
gas. The equation of state for such a system has the form: 

p = nT, n = ^ni. (18) 

Due to chemical freeze-out the conservation of particle number density currents for each species of hadrons is explicitly 



implemented in Eqs. ( |17[) . We ignore there the influence of possible changes of particle numbers due to decay of 
resonances to the time dependence of the macroscopic (hydrodynamic) variables, the validity of such an approximation 
will be discussed below in Sec. VI. 

We assume a quasi-inertial, Bjorken-like, longitudinal flow with velocity z/t dX the post-hadronization stage of 
evolution. Thereby 4- velocities are described by Eq. (Q). Using this flow ansatz and assuming that the hypersurfaces 
CT(a;), where the flow is in a quasi-inertial regime (^, correspond to the proper time T(r, 77) = \fW~^^ k, const, one 
can rewrite energy-momentum conservation equations ([l6| ) in the form ]4C| ] : 

v—]e+{e+p)D ^ 0, (19) 



dr &■ 



+ p) ( |: + «|:) + (|: + ) p - o, (20) 



where 



1 d 

-^. = 0, (21) 



D^l'\-+^-A^ + - + -- (23) 



dr dr 

In a similar way particle number conservation equations ( pT[ ) can be rewritten in the form 
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-^l^ + ^(7«)+7 - + - 



dr 



dr 



n, = 0. 



(24) 



Note that the equations (^j) - ( p4| ) are related, in general, to non- uniformly distributed energy and particle densities. 
Among them, the equation ( pi[ ) is the only one that governs the dependence of densities on fluid rapidity i]. The others 
contain 77 as a parameter. It allows one to analyze the solutions of the hydrodynamic equations ( p^ - (|2^) for any 
fixed rapidity. Our special interest concerns the point « corresponding to the main contribution to multiplicities 
at central rapidity. As we are interested primarily in the bulk properties of the system, we will search for the solutions 
of Eqs. ([1^) - ( p^ ) in the center of system, using a following approximations for transverse flow: 



V 

dv 
'd-r 



«1, 
dv 
dr 



(25) 
(26) 



The condition ( pSf ) is just a non-relativistic approximation for the transverse velocity (in LCMS for fluid element it 
coincides with v), whereas Eq. (|6|), as it will be discussed below, can be justified by a relatively large average energy 
per particle in a hadron resonance gas. Using the conditions (|2^ ) and (|2q), one can get from Eqs. (fL9|), (|2C| ) and (|2J) 
the approximate hydrodynamic equations 



d 

OT 

(e+p) 



d_ 

dr ^ 
d 
dr 
dv 1 
dr T 



P) 



dr 



dv 

dr 
d_ 
dr 

d 

d-r' 



= 0, 



p = 0, 



v-rr- I ni = 0. 
dr 



dr 



(27) 
(28) 
(29) 



Let us, motivated by the initial conditions at hadronization stage as well as exact analytical solution of non- 
relativistic hydrodynamics and Bjorken boost-invariant solution (see also ^]), assume the following ansatz 
for densities rii, transverse velocity v and temperature T in a central slice of space-time rapidity 77 sa 0: 



ni(r, r) = n^To, 0) 
p(t, r) = exp 

r 

v{t, r) 



^2 



p{T,r) 



T{T,r) 



Rv{t)' 
T{r). 



(30) 

(31) 

(32) 
(33) 



Here and below we omit indication of 77 - dependence of hydrodynamic variables. Using the above ansatz for the 
densities (|30|), ( ^l| ) and flow velocity distribution (^), we find that the particle number conservation equation (29) is 
satisfied if 



1 

Rv 



Rt 
Rt 



(34) 



Here and hereafter we designate by dotted symbols the derivatives in proper time: f— df /dr. 

Using the equation of state (|l8|) and conservation of particle number density currents (^) , one can rewrite Eq. (|2 
in the form: 



de{T) dT 

dT a7 



T 



dv 
dr 



0, 



where 



e(T) 



n 



e^iT) 



Hi 



'K2{m,/T) 



3T. 



Using the ansatz for transverse velocity distribution and temperature ( p2[ ) - (34), we get from Eq. 
ordinary differential equation for the transverse Gaussian radius parameter Rt and temperature T: 



(35) 
(36) 

a first order 
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Just for illustration one can solve Eq. ( ^ ) in the case of non-relativistic ideal Boltzmann gas when rui/T ^ 1. Then 
the energy per particle is equal to 



e,(r)«m, + 3/2r, (38) 



and one gets 



?2 ^__^__\ 2/3 



We will not use such an approximation which is too rough for real composition of hadron resonance gas. 

With the help of the equation of state (|T^ ) and the ansatz ( ^2|) - (34) for transverse velocity distribution and 
temperature the equation (28) is reduced to 

(e(T) + T) Rt Rt+ Rt RtT=t(^ + 2 - ^^^^ . (40) 

The last term in the bracket on the right hand side of the above equation is just — w^, as follows from Eqs. ( p2| ) and 
(^). Therefore, using the condition (p5|), one can neglect this term as compared with 1, and finally get the second 
order ordinary differential equation 

(e(r) + T) Rt Rt+ Rt RtT=T (^1 + 2 R^ +i^Z_^^ . (41) 

Together with Eq. (|3^) the above equation represents a system of ordinary differential equations for the transverse 
Gaussian radius parameter Rt and temperature T. These equations, together with ansatz ( ^^ - (^4|), correspond to 
the basic hydrodynamic equations ( p^ - (^4|) with ideal gas equation of state (O) and under conditions (p5|), (^6|) . 

In the above mentioned limit of the non-relativistic ideal Boltzmann gas {nii/T S> 1) the solution of Eq. { pT\ ) is 
given by Eq. ( ^ ) and so the (proper) time derivative of temperature is 

Then, for the same case of non-relativistic ideal Boltzmann gas, Eq. ( pl[ ) takes the form 

(43) 

If, in addition, one neglects the last two terms on the right hand side of the above equation , one arrives at 

T 

Rt Rt = ^ ■ (44) 

As is possible to show, Eqs. (|39| ) and ( ^4|) are exact equations for non-relativistic hydrodynamics with corresponding 
initial conditions (and obvious substitution t t). Here we just demonstrated the chain of steps which are necessary 
to pass from our semi-relativistic hydrodynamic equations ( |37|) and (^) for given initial conditions to non-relativistic 
ones. 

Let us now discuss the constraints for our ansatz parameters, that are forced by the conditions (p5|), (^). So far, 
as we will apply our solutions to the center part of system, limited by condition r/i?^ ^ 1, the first restriction ( p5| ) 
being applied to ansatz (^2|), ( ^4| ) means that 

i?T< 1- (45) 
The second constraint (p^), using ansatz (^2|), (Isl) and limitation r / Rt < 1, can be rewritten in the form 
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Rt Rj 



. 2 

Rt 



< 1. 



(46) 



and, taking into account 1^^), results in 
Finally, expressing Rt Rt by the help of Eqs. 



Rt Rt < 1 



and ([ll|), we get 



e(T) + T 



1 + 2rI+^^ 



(47) 



(48) 



We are interested in applying our solution to the post-hadronization stage of high energy heavy ion collisions, where 
typically r ^ Rt- Therefore, the above constraint is satisfied if the average energy per particle e{T) is large enough: 



T 



e(T) + T 



< 1. 



(49) 



The restrictions (^5|) and ( |49| ) define the region of applicability of our semi-relativistic solutions of hydrodynamic 



equations (19) - (24). 

Finally, we reduce the set of partial differential equations of chemically frozen hydrodynamics (^^ - (|2^) with 

(llSl) to the set of ordinary differential equations for the scale parameter Rt and 



relativistic ideal gas equation of state 

the temperature T , ( |37|) and (^), that can be solved by presently available numerical packages, e.g., by the program 
package Mathematical^ . The initial conditions for these equations are: Bjorken-like initial longitudinal flow, where 
velocity z/t is accompanied by transverse flow, and Gaussian density profile in transverse direction. Such reduction 
is based on approximations (psl), (|26|), which lead to the semi-relativistic hydrodynamic equations (|27| ) - (^9|), and, 
as we will see, are satisfied in central domain r / Rt 1 at post-hadronization stage of matter evolution in Pb +Pb 
collisions at SPS CERN. 

One can see that the equations (|3^) , ( pl| ) depend only on the concentrations Ki of different particle species and do 
not depend on absolute values of densities, and, thus, do not depend on the unknown common chemical potential fXch- 
The concentrations k.^ are completely defined by the temperature and baryochemical potential that are extracted from 
the fitting of particle number ratios (see Sec. III). Then Eqs. (^7|), (plj ) give the possibility to find the proper time 
of chemical freeze-out Tch if one knows the "initial conditions" at thermal freeze-out Tth'- temperature T, transverse 
scale parameter Rt and velocity at the Gaussian border Rt- These values can be found from analysis of pr-spectra 
and HBT correlations. Then solving (^7|), ( pl| ) numerically with initial conditions at thermal freeze-out, we find T{t), 
and then find Tch from the equation Tjrch) ~ Tch- It gives the possibility to find the effective volume (|l^) at chemical 
freeze-out, V^j^j- Further, using the absolute value of the rapidity density for some particles, e.g., for pions, we can find 
the common chemical potential fich and, so, the common fugacity - the factor that changes the energy and particle 
number densities as compared to chemically equilibrated ideal gas. Finally all densities at chemical freeze-out can be 
determined. Using the Eqs. (p2), (p4) one can also obtain the transverse velocity u(r, r) at chemical freeze-out. 



V. RECONSTRUCTION OF THERMAL FREEZE-OUT STAGE 



To reconstruct the thermal freeze-out conditions, we assume that the freeze-out happens in a rather narrow time 
interval at approximately constant temperature Tth- It corresponds in our solution to some proper time Tth ~ const - 
Of course, the different particle species can have different freeze-out (proper) times, so the time Tth is just some mean 
value, when majority of strong interacting particles (pions, protons, etc.) leave the system. Such an approximation 
gives no big influence on the reconstruction of hadronization stage because, as we will show, the deviations of different 
freeze-out times from that mean value for majority of particles are relatively small. As for hyperons, like fJ, that 
can escape from the system essentially earlier, it also does influent a little on hydrodynamic evolution because of a 
relatively small numbers of those particles. 

Within our fltting procedure we do not use temporal widths of thermal freeze-out hypersurface because of two 
reasons. 

First, such a width is introduced usually just as a common pure temporal factor to the local equilibrium functions. 
It contradicts to natural picture of emission, where the space-time correlations have to be presented: particles from 
a periphery escape more easily than from the center (see discussion about it in Ref. p^). The second point is that 
such a temporal smoothing of locally equilibrium emission function is incorrect, at least, if one based on Boltzmann 
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equation (BE) For example, in order to calculate spectra in the case of exact locally equilibrium solution of BE, 
one can use either Cooper-Frye prescription ||2^ (sudden freeze-out) or, alternatively and equivalently, continuous 
emission with non- equilibrium form of emission function as it is shown in Ref. pq] . 

We choose the transverse profile factor in Gaussian form, p(r, r) — exp{—r^/2R^{T)), as we noted before. Then we 
have for the effective volume (|TT|): 

Veff « 27TTRUr). (50) 

To extract the temperature and transverse flows at thermal freeze-out we fit in our model the experimental spectral 
slopes for different particle species. 

The transverse spectra for direct particles can be expressed through the Wigner function: 

- f da^i^W! { '-4^, ^) Pi^y (51) 



2'KmTdmTdy J ^ \ Tth ' Tth 

At a region of small px the contribution from the decay of resonances has a peak that results in a significant change 
of the slope of transverse spectra. At large pt , the relative contributions of a resonance decays are distributed much 
more uniformly in pT, especially if transverse flow takes place, so, while spectra are shifted to higher values due to 
resonance decays, the spectra slopes do not change essentially and the slopes of the spectra are approximately the 
same as for direct particles ■ Below we neglect the contribution of resonance decays into the spectra slopes in that 
region. At very high energies pr > 2 GeV particles from the hard partonic processes at the early collision stage may 
contribute to the spectra. We are fitting the transverse spectra according to Eq. (|5l| ) in the mr— interval between 
0.3 and 2 GeV/c for pions, kaons and protons. For hyperons we use an interval up to 3 GeV/c. 

The other important point is how our results are influenced by the velocity profile. The linearity of the transverse 
velocity in r ( ^2| ) is associated with the non-relativistic approach for transverse flow ( p5[ ) and is valid for r < Rt ■ 
Note that the integration in Eq. ( pT| ) is carried out over all r-region requiring to use correct velocity distributions 
when u(t, r) < 1 in the whole region r < oo. So we have to extrapolate relativistically the velocity distribution in 
the r-region of applicability of our hydrodynamic solution, where w^(t, r) 1, to larger r . For this aim we use the 
results of Refs. |^,^ where two "extreme" relativistic extrapolations of a linear transverse velocity distribution from 
non-relativistic region to relativistic one was used: 

viTth,r) ^tanhirv'{Tth,0)) (52) 

and 



v{Tth,r) = rv'{Tth,0)/V^ + irv'{rth..OW, (53) 

where v{Tth, r) « rv' [rth, 0) for small transverse velocity v{Tth, r). The extrapolation ( p2[ ) corresponds to "hard" flow, 
while the extrapolation ( ^3| ) - to "soft" one |^ , ^ . We have used both of them to fit spectra and found that the 
deviation in extracted temperature Tth as well as in v' [tih, 0) is about 5 percent. It means also that for SPS energies the 
relativistic tails of transverse velocity profile have no essential influence on the spectra in above mentioned momentum 
region, the transverse expansion, therefore, is basically non-relativistic. Thus we reproduce the plots of spectra for 



only one velocity profile corresponding to the "hard" flow, see also discussion in Refs. 43 45|]. 

There are two fitting parameters in the analysis of the slopes of transverse spectra according to Eq. ( pT| ) : Tth and 
o-th = {v'{Tth,0)RT)~^ , where is the intensity of transverse flow p3| , ^ , p7t 

^y— {v'rp{T,0))^^ hydrodynamical length ^^^^ 
Rt transverse radius 

One can see that 1 / ^/a is the hydrodynamic velocity at the Gaussian " boundary" of the system, 

VB^riT) = '"{t^''^ ^ Rt{t)) = (55) 



within a linear approximation for the velocity profile, and according to Eqs. ([32|), (34) 



1/V^= Rt. (56) 

For very intense relativistic transverse flow, a — > 0. On the other hand for a 3> 1 we have non-relativistic transverse 
flow. The results of the fit are presented at Fig. 1 and Fig. 2. 
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FIG. 1. The one-particle transverse mass spectrum for p — p h~ and Kg Solid line corresponds to Tth = 135 MeV 
and ath = 7, dashed line: Tth ~ 100 MeV and Uth ~ 4. The overall normalization is arbitrary. 
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FIG. 2. The one-particle transverse mass spectrum for D, = D,^ + Q , H~ and A Solid line: Tth = 135 MeV and 

ath = 7, dashed line: Tth = 100 MeV and Uth ~ 4, dash-dotted line correspond to chemical freeze-out parameters: Tch = 163.37 
MeV and Uch = 10.8. The overall normalization is arbitrary. 

As one can see from Fig. 1 the pion spectra can be described with good accuracy in two regimes: relatively low 
temperature of thermal freeze-out, Tth = 100 MeV and relatively intensive transverse flow, ath = 4, on the one hand 
and higher temperature Tth — 135 MeV and non-relativistic transverse flow (within the typical size of the system 
at freeze-out Rt) , ath ~ 7, on the other hand. While the latter parameterization describes well the spectra of all 
particles, except, possible, the baryon 17, the former cannot reproduce the data on spectral slopes except for pions 
(see Figs. 1, 2). 

It is worth to note that the transverse spectra for Q, could be described rather well, if one assumes that it escapes 
from the system just after hadronization with corresponding chemical freeze-out parameters: Tch — 163.37 MeV 
and ach — 10.8 , the latter is evaluated from the hydrodynamic equations (see next section) with "initial" values 
Tth — 135 MeV and att = 7 at thermal freeze-out. The satisfactory result for the O -transverse spectra with chemical 
freeze-out parameters, presented in Fig. 2, justifies the assumption of earlier freeze-out for f2 that could happen just 
after hadronization due to small cross-sections of these particles. 

Finally we fix the thermal freeze-out parameters Tth = 135 MeV and ath = 7 according to the best fit for the 
set of transverse spectral slopes for pions, kaons, protons. A, S and 51. The reason to use another, low tempera- 
ture parameterization |^, was that relatively large transverse flow gives the possibility to describe satisfactorily the 



decrease of the observed pion interferometry radii with increasing of 18 |. As a consequence, the temperature of 
thermal freeze-out becomes low, about 100 MeV. Such interplay between flow and temperature is easy to understand 
qualitatively from a simple approximation for spectral slopes: 



Ti,eff ~ Tth + mil ath, 



(57) 



^Note here, to avoid the misunderstanding, that the thermal freeze-out parameters T — 100 MeV and a = 4 do not correspond 
to the our solution of hydrodynamic equations with the initial conditions at T = 135 MeV and a = 7. 
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that is valid for heavy mass particles and non-relativistic flow [|7j. However, as one can see from Fig. 1 and Fig. 
2, such a low temperature parameterization fails to reproduce the spectra slopes simultaneously for many particle 
species, as we demonstrate for p, K, A, S and fl particles. 

To discuss the problem, we would like to draw attention to resonance decay as an important factor that could 
influence the behaviour of the interferometry radii in px , sl factor that was not yet taken into account quantitatively. 
It is rather difficult because one needs to take into account of a large number of resonances. We note that the wide 
spread opinion that due to relatively small temperatures in hadron resonance gas one can use only resonances with 
low masses is not quite correct and, at least, needs verification, because the number of resonance states at a fixed 
mass interval increase quickly with mass. 

Generally, the contribution of resonance decays into the final particle correlation function is a rather complicated 
effect. Namely, the resulting resonance contribution to the interferometry radii depends mostly on the following: 
statistical weight and mass distribution of resonances in the rest frame of fluid element, the intensity of flow, lifetime 
of resonance, feeding from middle-lived resonances, such as w, and the kinematic "window" for pions coming from 
resonance decay. The important factor is also the length of the region of homogeneity , that is effective size of 
fluid element emitting the resonance in some px region, for developed flows it is rather less than the corresponding 
values for pions due heavier resonance masses. All these factors vary for different resonances and we cannot a priory 
exclude or easily estimate the influence of resonance decays on the interferometry radii at CERN SPS A+A collisions. 
Qualitatively, the steep decrease of the observed interferometry radii on pr could be due to an interplay between two 
factors: transverse expansion of the hydrodynamic tube and contribution to the radii from resonance decays. For 
relativistic transverse flow the contribution of low-mass resonance decays to the interferometry radii is negligible [Q , 
iQ and the decrease is determined mostly by flow. If the transverse expansion is not very intensive, the decrease 
of the interferometry radii on px is the result of both effects: flow and decays of resonances |Q. Therefore, in 
our opinion, the detailed evaluation of the resonance contribution to interferometry radii in hydrodynamic approach 
needs a serious study, taking into account the contribution of high mass resonances. In this paper we analyze the 
interferometry radii for pions, kaons and protons at maximal measured pr region only, in order to minimize the 
influence of resonance contribution. Here we are based on known results |^ , ^ of hydrodynamic approach and our 
preliminary studies with relatively small number of resonances ji^ . According to these results, the resonance decays 
lead to an increase of the interferometry radii at small pT < ptq ( Pro ~ 0.3 0.4 GeV for pions) and to a vanishing 
of resonance contributions into the radii out of this region. The latter means that the complete interferometry radii 
and interferometry radii for "direct" particles (calculated neglecting resonance decays) are rather close to each other 
at pt > Pro- 

Finally, we would like to note that the aim of the paper is not the description of all spectra and correlations in 
the whole momentum region but the reconstruction of hadronization conditions. For this aim it is enough to use 
the slopes of spectra only and interferometry radii in some momentum interval. As for the description of spectra 
and correlations in whole momentum region, it is a separate problem that demands a detailed analysis within our 
hydrodynamic model of decays of huge number of resonances, and taking into account the back reactions to these 
decays at all stages of evolution of hadronic system (to calculate the number of resonances which are survived till the 
thermal freeze-out). 

Our calculation of Bose-Einstein and Fermi-Dirac correlation function, 



is based on the distribution function f {Pfj,uf^ (r, rj) /Ttii{ri) , fii tii(r]) /Tchiri)) in Boltzmann approximation at the hyper- 
surface of thermal freeze-out. For tt^tt" correlation function we have, e.g., 







th ' '^th J 


/9(a:)e*?^ 


2 


[jdafI^{x)p^j{E^,^) 




(/d<(-K/(2^,%f) 





c{p. 9) = 1 + A , ... \ " „ ":^ .. ... . (59) 



where p = (pi +p2)/2, q = pi —p2 and A is phcnomenological parameter describing the experimental suppression of the 
correlation function due to, e.g., contribution of long-lived resonances. We are fitting by -method the correlation 
function with fixed parameters T^^ = 135 MeV and ath = 7 to Gaussian approximation of experimental data: 
C{p, Qi, 0,0) = 1-1- Xexp[—qfRlf{p)] in Bertsch-Pratt parameterization. Then i = I, s, o correspond to longitudinal 
(I) and two transverse (t): sideward (s) and outward (o) directions. 

Our results for proper time Tth and geometrical Gaussian radius Rt of thermal freeze-out are presented in Table 3 
in the last two columns. The experimental interferometry radii Ri and Rs for pion, kaon and proton pairs measured 
in LCMS near mid-rapidity are taken at maximum possible px ■ We assume, in agreement with experimental data for 
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pion distribution in rapidity, that the factor pir, r, rf) does not vary strongly in rapidity within one unit of rapidity 
near rj = Q. 

TABLE 3. The characteristics of the thermal freeze-out in Pb+Pb SPS collisions. ^ 





Pt GcV 


Ri fm 


Rs fm 


R.X fm 


Tth fm 


TT TT [50| 


0.5 


4.6 


4.25 


5.3 


8.2 


K+ K+ 151] 


0.91 


3.2 


3.59 


5.1 


8.9 


PP |52| 


1.073 


2.9 


3.4 


5.2 


9.5 



Let us note that the proper time of thermal freeze-out Tth and transverse radius Rt of system is slightly different 
for different pair correlation functions. To reconstruct the chemical freeze-out we will use the average values of the 
above parameters. In the Table 4 we compare our results (1st fit) with the results of other authors. Temperature Tth, 
proper time Tth, transverse Gaussian radius of the system RxiTth) a-ud transverse velocity at Gaussian "boundary" 
of the system, Vfij,{Tth)— ^/^/cuth, at thermal freeze-out for fits 2 — 4 are taken from Refs. ||l8|, and for fit 5 from 
Ref. [ p^ . Note, that in our model of post-hadronization hydrodynamic evolution the transverse velocity at Gaussian 
"boundary" of the system coincides with proper time derivative of transverse Gaussian radius, see Eq. (^6|). In the 
columns 4, 5 we present our results for the values of chemical potential of direct pions fj.^^th and particle number 
density (before resonance decay) nth = n{Tth,r = 0, 77 = 0) at the thermal freeze-out in the very center of the system, 
which were calculated by the method described in the next section for different parameters Tth, nh, Rt {Tth), VR^{Tth) 
(fits 1-5). 



TABLE 4. The characteristics of the thermal freeze-out in Pb+Pb SPS collisions. 



Fit 


Tth MeV 


Tth fm 


RriTth) fm 


VRriTth) 


nth 1/fm^ 


Mtt,*/! MeV 


1 


135 


8.9 


5.2 


0.378 


0.274 


74.3 


2 


100 


7.8 


6.5 


0.424 


0.199 


128.5 


3 


120 


6.6 


5.9 


0.339 


0.286 


119.7 


4 


160 


5.5 


5.6 


0.247 


0.381 


50.0 


5 


139 


5.9 


7.1 


0.55 


0.221 


32.9 



At the first sight, the particle number density at the thermal freeze-out nth = 0.274 1/fm^, which is obtained in our 
analysis, is too high, formally, the mean free path for pions, as we will demonstrate below, is rather small as compare 
with transverse "size" of the system, RriTth)- For static systems such a result could indicate the absence of the 
freeze-out of momentum spectra at this temperature. However, if the system expands intensively, the comparison of a 
rate of expansion with the rate of collisions plays the dominant role for freeze-out criterium 1 53 - 56 1 , , and is more 
important for estimation of freeze-out parameters than the comparison of mean free path with geometrical Gaussian 
"size" of the system. Namely, for the system that we are interested in, local thermal equilibrium and ideal fluid 
hydrodynamic behaviour can be maintained as long as the collision rate among the particles is much faster than rate 
of expansion. Since densities drop out during 3-dim expansion rather intensively, the collision rate decreases rapidly 
and at some moment system falls out of equilibrium. As a result, system decouples and the freeze-out happens. Then, 
to describe momentum spectra we use the standard Cooper- Frye freeze-out prescription [2^ ] , that assumes validity of 
the ideal hydrodynamics up to a "sharp" 3-dim freeze-out hypersurface, whereas outside of this surface the particles 
are supposed to move freely. It is worth to note, that Cooper-Frye prescription can be applied only at densities where 
the matter is (locally) in thermal equilibrium, and therefore at the thermal freeze-out the rate of expansion should be 
comparable with the rate of collisions but smaller than it. Let us check, whether the densities, that correspond to our 
fitting parameters. Table 4 (fit 1), meet to such a freeze-out criterium. 

The inverse expansion rate is the collective expansion time scale |5q]: 



1 dn 
n dt* 



u'^df^n 



1 



(60) 



where t* is proper time in the local rest system of fluid and the last equality follows from the conservation of particle 
number density currents ([l7|). Using for 4- velocities Eq. (^, one can easily get that 



^Two-dimensional fits for kaon and proton correlation functions in and were used . 
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2i?T 



n ,, d d , . ,1 , ^ 

5m"" = ^7+^(7^) + 7(- + -). (61) 

Taking into account expression for transverse velocity (|3|), (|3^), we get finally 

(62) 

Using our parameters at the thermal freeze-out from Table 4 (fit 1) we obtain that collective expansion time at thermal 
freeze-out Tc^pinh) « 3.88 fm. Note, that Texp{Tth) < Rrinh) = 5-2 fm. 

The inverse scattering rate of particle species i is the mean time between scattering events for particle i, t^.^^^.: 

(63) 




where Vij is the relative velocity between the scattering particles and cTy is the total cross section between particles 
i and j, the sharp brackets mean an average over the local thermal distributions. This time is determined by the 
densities of all particles with which particle i can scatter, and the corresponding scattering cross sections. Let us 

estimate the mean time between scatterings for pions, t^^I^. First, note that r^^^^jj > A'^*^ where A^*^ is mean free path 
for particle species i: 

A« « , (64) 

so Ag*2i( represents the lower limit for t^^J^^. Let us estimate the pion mean free path in the very center of the system, 
^'"^ {Tth, r = 0), using the following approximation: 

X^^\m.,0) « i (65) 

where CTttp = 65 mb is supposed to be total cross section for pion-proton scattering, and we assume the same cross- 
section for all non-strange baryons, whereas aT^-n = 10 mb is total cross section for pion-pion scattering and we 
assume the same cross-section for all non-strange mesons. The densities ub = 0.044 l/fm'^ and um = 0.159 1/fm^ 
are the densities of non-strange baryons and non-strange mesons respectively at r = Tth and r = 0, the densities 
were calculated by the method described in the next section. The contribution of strange particles to pion mean free 
path is neglected in (|65|). Then X^'^\Tth,0) ~ 2.26 fm and this value is essentially smaller than RriTth) = 5.2 fm, 
while the mean time between scattering events for pions rfj^t > X^^^ and, therefore, is less but comparable with the 
collective expansion time Tcxp{Tth) ~ 3.88 fm. So, one can conclude that our fitting parameters correspond to freeze- 
out criterium. Note, that the pion mean free path at the Gaussian "boundary" of the system A^'^^(rt/i, Rt) ~ 3.73 fm, 
and, thus, the densities, that are represented in Table 4 (fit 1), are just at the border where ideal hydrodynamics still 
can be applied. Lower values for densities, on our opinion, are not appropriate to apply the Cooper- Frye prescription 
within hydrodynamics. 

Of special interest is the average transverse velocity, that depends on the equation of state during the thermalized 
matter evolution after collision. At central slice of hydrodynamic tube the mean transverse hydrodynamic velocity 
{v) = {\vt{ti — 0)1) and mean square of transverse hydrodynamic velocity (w^) are defined as follows: 



^ /$l^/^P^^"(^)/(x,p) ^ j'-^u^{x)p{x)v-{x) 
In the first approximation in 1/a expansion one can get for a ^ 1. 



1,2. (66) 



f, (67) 
2a a 

This approximation corresponds to non-relativistic flow and is model-independent (the result does not depend on 
relativistic extrapolation of the velocity profile, see discussion above). Then, at thermal freeze-out {ath = 7) we get 
for the average transverse velocity of the central (j] = 0) slice of the hydrodynamic tube: (u)^^ « 0.47. 

Note, that using (E^) , one can get for cylindrically symmetric expansion the expression for transverse spectral slopes 
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1 1 2 

Tt,eff ~ Tth + -rui {vlf^) = Tth + —rrii (vth) ■ (68) 

Right-hand side in the first equahty is approximately the ratio of average transverse energy of particles species " i" in 
central slice, ^ ^"^^^ ^ ~ ^^Tth + ^^"^-i {'"th) to total number of particles of such species in central slice, at 
thermal freeze-out. 

VI. HYDRODYNAMIC RECONSTRUCTION OF CHEMICAL FREEZE-OUT STAGE 

Based on the hydrodynamic solutions found in Sec. IV and "initial conditions" at thermal freeze-out considered in 
Sec. V, we can restore the space-time extension of hadronization stage and find the particle and energy densities. We 
use the solution of hydrodynamic equation in the center of the system, where the majority of particles are produced and 
where non-relativistic approximation for transverse velocity is correct. In Fig. 3 we demonstrate the cooling down of 
HG as it depends on proper time r for fit 1. The "initial condition" for the hydrodynamic solution is the temperature 
of thermal freeze-out Tth = 135 MeV at corresponding time Tth — 8.9 fm/c which were found from the spectra 
and correlations in the previous section. The "initial" velocity ( at Tth) at the Gaussian "boundary" of the system 
VB.T {'''th) ~ 0.38 justifies the applicability of our approach, namely, non-relativistic transverse expansion ( p5|) in central 

domain: v\^{Tth) =Rt (Tth) ~ 0.14 <C 1 and so the restriction (^ ) is also satisfied between chemical and thermal 
freeze-out. The second condition ( p^ ) and restriction (|4^), that follows from it, is also satisfied: Tch/ {e{Tch) + Tch) ~ 
0.141 < 1 and Tth/{e{Tth) + Tth) ~ 0.127 < 1. 
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0.12 ^1 1 ^ 1 ^ 1 ^ 1 ^ 1 1 
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T (fm) 

FIG. 3. The temperature T (GeV) of the hadron resonance gas as function of the proper time r (fm/c) of its evolution. 
Sohd line corresponds to the semi-relativistic solution T{t) of Eqs. ( ^7| ) and (^). Dashed line describes T{t) in completely 
non-relativistic picture of evolution of Boltzmann gases, Eqs. ( p^ and (|44[), with the same initial conditions. The effective 
mass of the mixture of non-relativistic Boltzmann gases is m = ^ Kirrii « 0.662 GeV. 



One can restore the proper time of the chemical freeze-out from the intercept of the curve T(t) at the temperature 
of the chemical freeze-out T^h ~ 163.37 MeV, that was found from the analysis of particle number ratios. Using the 
result, Tch — 7.24 fni/c, one can calculate the effective volume V^j:''j ~ 2TiTchR^{Tch) ~ 975 fni^ associated to unit 
of rapidity at mid-rapidity, and then find common chemical potential ^Lch from pion multiplicities at mid-rapidity, 
dN^- /dy « 160 For direct pions the common chemical potential fj,ch is added to phenomenological one: /ijr.ch = 
Mtt + ^J'ch■ Then one can easily evaluate the energy and particle densities. 

Before proceeding to numerical results, let us check the validity of our basic approximations. The duration of 
hydrodynamic expansion of hadron gas. At — 8.9 — 7.24 = 1.66 fm/c, is smaller or comparable with the relaxation time 
to chemical equilibrium in the system (see, e.g., [ |57| , p8t ), that exclude the chemical equilibration during the evolution. 
The typical life-time of the short lived resonances, is, however, compared with the period of hadron hydrodynamic 
stage At — 1.66 fm/c, therefore, such decays could take place at the background of chemically non-equilibrium 
expansion. The decays could result in an increase of number of (quasi) stable particles during the expansion (without 
that all particle numbers will conserve within concept of chemical freeze-out). On the other hand, there are back 
reactions, that produce the resonances. These reactions do not lead to chemical equilibrium, but just reduce effect 
of resonance decays into the chemically non-equilibrium medium, resulting in an increase of effective life-time of 
resonances. 

At the same time, some residual effect of resonance decays during the evolution could, anyway, take place. We ignore 
that effect since our approximate calculation of resonance decays demonstrate two opposite influences of it on a fluid 
evolution. (To do this estimate we replace the discrete spectra of heavy resonances by the Hagedorn distribution). 



18 



The decay of heavier resonances with subsequent thermalization of decay products leads to a heating of the system, 
while the decay of relatively light resonances results in cooling. The reason is clear: the difference between kinetic 
energy of decay products and kinetic energy of thermal motion of the same particles (e.g. pions) is positive in the 
former case and negative in the latter one. Hence, the corresponding effects are approximately compensated and give 
no big influence on hydrodynamic evolution. 

Our results are presented in Table 5. The first line corresponds to parameters at thermal freeze-out which are 
represented in the first line of Table 4. In Table 5 we present also our results for proper time and densities at the 
hadronization stage based on the thermal freeze-out parameters taken from other papers: fits 2-5 from Table 4. 

From the results presented in Table 5 one can see that the ratio of energy density to total particle density at 
chemical freeze-out reaches the value e chinch = 0.997 GeV, in agreement with the phenomenological observation in 
Ref. 1^. The energy density at the hadronization point 'EcK ~ 0.42 GeV/fm"^ is close to results of lattice calculation at 
the region of the phase transition and about equal to the energy density inside a nucleon. Note that for temperature 
and chemical potential found in |^ at the chemical freeze-out, we get Sch ~ 0.54 GeV/fm'^ and Sch/nch ~ 1-16 GeV, 
that is close to our results. Using the values of concentrations of HG and temperatures at chemical and thermal 
freeze-out, we get for conventional velocity of sound Cq = p/e = T/ J2 Kiei{T) the values 0.164 and 0.145 at chemical 
freeze-out and thermal freeze-out, respectively. Such relatively small values are determined by the large average mass 
of particles in hadron resonance gas, that is 0.662 GeV for hadronic masses below 2 GeV. 

In Table 4 we list also the total particle number density nth = richV^jj /V^^j, and chemical potential for direct 
pions fJ,Tr,th at thermal freeze-out. One can see that for small temperature fits, Tth = 100 120 MeV (fits 2 and 3) 
the chemical potential of direct pions at the stage of system decay seems to be too large and closed to condensation 
point, fi^,th = m^- 

Using Eqs. (|56|), ( |67| ) and velocity at Gaussian "boundary" of the system at the chemical freeze-out f_Rr ('''c/i)=i?T 
{Tch), we calculate the averaged velocity at chemical freeze-out (f)^^ for fit 1 and find that (w)^^ ~ 0.383. 



TABLE 5. The characteristics of the chemical freeze-out in Pb+Pb SPS collisions. 



Fit 


Tch fm 


fich MeV 


RriTch) fm 


VRT{Tch)=RT {Tch) 


Ech GeV/fm3 


Hch l/fm^ 


1 


7.24 


-67 


4.63 


0.304 


0.420 


0.421 


2 


4.1 


-9.8 


5.18 


0.278 


0.598 


0.600 


3 


4.2 


-13.5 


5.19 


0.245 


0.583 


0.585 


4 


5.3 


-76.1 


5.55 


0.240 


0.398 


0.399 


5 


4.8 


-110.4 


6.53 


0.491 


0.323 


0.324 



VII. CONCLUSIONS 

The reconstruction of the hadronization stage for CERN SPS Pb-|-Pb collisions at 158A GeV/c is obtained within the 
hydrodynamic approach for the post-hadronization stage. Because of incomplete stopping at SPS energy we demon- 
strate that the analysis of particle ratios can be done within a relatively narrow rapidity window near mid-rapidity, 
provided one takes a non-uniformity in mto account. We estimate that the non- uniformity of baryochemical 
potential may result in small but noticeable (w 10 %) contribution to the pion multiplicity at mid-rapidity due to 
transport of pions from resonance decays from non-central rapidity regions. Using an optimization procedure for the 
calculation of particle number ratios near mid-rapidity, and also reconstructing the thermal freeze-out stage from 
spectra and correlations for different particle species, we find the physical conditions at the hadronization stage. The 
proper time for chemical freeze-out is found to be too small, Tch ~ 7.2 fm/c to reach the chemical equilibrium in pure 
hadron gas, or cascade "event generator" , approach. We stress that our method does not assume complete chemical 
equilibrium of Boltzmann gases at the hadronization, and allow to estimate the deviation from it. We find that 
distribution functions of hadrons at hadronization hypersurface include some common factor, exTp{fj,ch/Tch) ~ 0.66, 
that reduces essentially particle numbers as compare to chemically equilibrium ones. The energy density is found to 
be close to the results of lattice calculations at the point of the phase transition. The ratio of energy density to total 
particle number density is found to be about 1 GeV. 

All those results indicate the existence of a pre-hadronic stage of expansion before the chemical freeze-out and the 
comparison with lattice calculations pointed to a QGP phase as most probable state of the matter till r ~ 7.2 fm/c 
in Pb+Pb 158 A GeV/c collisions. The obtained characteristics of the hadronization give the possibility to make the 
next step and reconstruct the pre-hadronic stage of evolution and, then, to study the basic properties of QGP. 
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